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Abstract 

A dressing of a nonspherical potential, which includes n zero range potentials, is 
considered. The dressing technique is used to improve ZRP model. Concepts of the 
partial waves and partial phases for non-spherical potential are used in order to perform 
Darboux transformation. The problem of scattering on the regular X n and YX n structures 
is studied. The possibilities of dressed ZRP are illustrated by model calculation of the 
low-energy electron-Silane (S1H4) scattering. The results are discussed. 
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1 Introduction 

The ideas of zero range potential (ZRP) approach were recently developed to widen limits of 
the traditional treatment by Demkov and Ostrovsky QXJ] and Albeverio et al. [0. The advantage 
of the theory is the possibility of obtaining an exact solution of scattering problem. The ZRP 
is conventionally represented as the boundary condition on the matrix wavefunction at some 
point. Alternatively, the ZRP can be represented as pseudopotential (Breit ||). 

On the other hand, Darboux transformation (DT) allows to construct in natural way exactly 
solvable potentials. General starting point of the theory goes up to Matveev theorem (see [|J). 
The transformation can be also defined on the base of covariance property of the Schrodinger 
equation with respect to a transformation of wavefunction and potential (Matveev and Salle 
||). Darboux formulas in multi-dimensional space could be applied in the sense of Andrianov, 
Borisov and Ioffe ideas 0. In the circumstances, DT technique can be used so as to correct 
ZRP model. 



We attempt to dress the ZRP in order to improve the possibilities of the ZRP model. We 
use notations and some results from ||13|| . DT modifies the generalized ZRP (GZRP) boundary 
condition (Section |2|) and creates a potential with arbitrarily disposed discrete spectrum levels 
for any angular momentum I. In the section |^ we consider S'-representation for a non-spherical 
potential so as to dress a multi-centered potential, which includes n ZRPs. As an important 
example, we consider electron scattering by the X n and YX n structures within the framework of 
the ZRP model (Section f|). In section [O] we present the our calculations for the electron-SiH 4 
scattering and discuss them. 



2 Zero range potential dressing 



Let us start from the simplest case of a central field. Then angular momentum operator com- 
mutates with Hamiltonian and therefore wavefunction \I/(r) can be expanded in the spherical 

waves 

oo I 

tt(r) = Mkr)" 1 J2 E e*to +l "M ^(r, fc)Y, m (n)Y^(n Q ), (2.1) 

1=0 m=-l 

where n = r _1 r, n is initial particle direction, ipi(r,k) are partial waves, and rji are phase 
shifts. Consider the radial Schrodinger equation for partial wave with angular momentum I. 
The atomic units are used throughout the present paper, i.e. h = m = 1 and Born radius 
ao = h 2 /me 2 = 1. 

H^i(r,k) = k 2 ^(r,k), (2.2) 

H l = -D 2 + 2u l (r), (2.3) 

D = d/dr denotes differential operator, and Hi are Hamiltonian operators of the partial waves. 
This equations describe scattering of a particle with energy E = k 2 /2. The wavefunctions 
ipi(r, k) at infinity have the form 

ipi(r, k) r ~°° sin \kr - y + j . (2.4) 

Let us consider GZRP in coordinate origin. This potential is conventionally represented as 
boundary condition on the wavefunction (see [0]) 



(r%(r,k)) = - , , ^, (r%(r,kj) , (2.5) 

V r(V ' >)r=o (21 + 1)!!(2Z - 1)!! ^ J 'r=o ' v ' 

where ai are inverse scattering lengths. The potential u\ = 1(1 + l)/2r 2 and therefore wavefunc- 
tions ipi(r, k) can be expressed in terms of the spherical functions 

if)i(r, k) = cos(r]i)ji(kr) + sm(rji)ni(kr) , (2.6) 

where spherical functions ji(x),ni(x) are related to usual Bessel functions as ji(x) = y^f Jz+i/2(^) 
rii(x) = (—\) l ^p^-J^i^i/2(x). In the vicinity of zero they have the asymptotic behavior ji(x) ~ 

x l+l /(2l + 1)!!, and rti(x) ~ x~ l (2l — 1)!!. To substituting the equation (|2.6| ) into the boundary 
condition we obtain the elements of S'-matrix 

aci — ik 2l+1 



2 



The bound states correspond to the poles of the S'-matrix (i.e the zeros of the denominator 
ai + ik 2l+1 ), which lie on the imaginary positive semi-axis of the complex fc-plane. It is obvious 
that bound state, with orbital momentum I, exists only if (— l) z otj > (elsewise an antibound 
state exists) and has the energy E\ = a 2 ^ 2l+1 ' /2. 

Thus, spectral problem for GZRP is solved for any value k. On the other hand, the equations 
( ]2.2| ) are covariant with respect to DT that yields the following transformations of the potentials 
(coefficients of the operator Hi) 

ui 1] =ui- Ds h (2.8) 

and the wavefunctions ipi(r, k) 

W'> = ^=£=AM), s l = ^^ (2.9) 
<k 2 + b? <Pi{r,ibi) 



where ipi(r,ibi) are some solutions of the equations ( [2.2[ ) at k — ibi, and bi are real parameters, 
which can be both positive or negative. The DT (|2.9| ) combines the solutions ipi{ r ,k) and a 
solution ifi(r,ibi) that corresponds to another eigen value k 2 = —bf. Repeating the procedure 
we obtain a chain of the integrable potentials u\ . In general, dressed potential u\ is real for 
real function ipi(r,ibi). 

The next step in the dressing procedure of the zero-range potential (ui = 1(1 + l)/2r 2 ) is a 
definition of the free parameters of the solutions <fi(r, ibi). Suppose the prop functions ipi(r, ibi) 
satisfy the boundary conditions (|2.5[) with ai = e\. In the simplest case of I = we have 



(fo(r,ib) = sinh(6r) — ib/e) cosh(6r), (2-10) 

and 

So ~ — e + [b — e Jr. 
The DT ( |2.9| ) gives rise to the following requirement on dressed wavefunction 

\ a — e I 



The dressed potential is given by 



„(i) 



b 2 {b 2 -e 



2\ 



(6cosh(6r) — esinh(6r)) s 



It is regular on semiaxis only if (6/e)e[0, 1). In the limiting case at b —>■ we obtain long- 
range interaction oc (1 — er)~ 2 , which can be regular on semiaxis only if e < 0. Assuming 
e = ±6 we get = (trivial transformation), and boundary condition can be obtained by 
the substitution: 

b 2 + k 2 ^ h 

a — > — ± o. 

a =F b 

To dress free wave (a = oo) we obtain ZRP at the coordinate origin. Thus, ZRP can be also 
introduced in terms of DT. To consider transformation with parameter e = a we obtain regular 
solution ifj^ and tangent of phase shift is 

tanr ?o " be + k 2- 
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In the other cases asymptotic of the functions si at zero is given by 



r-»0 , 



I 



b 2 r 



r (21-1)' 

l + l 
+ 



t\ 7^ oo, I > 



b 2 r 



(21 + 3) 



ei = oo. 



It is clear that the each DT introduces short-range core of centrifugal type (which depends 
on angular momentum I) in the potential. In this situation the boundary conditions on the 
dressed wavefunctions require modification. Thus, in the case e\ ^ oo, / > the boundary 
conditions become 



r=0 



and in the case e\ = oo we obtain 



or,(2Z-l)! 



(k 2 - b 2 )(2l + l)\\(2l - 



r=0 



D 



V Vl )r=o (k 2 -tf)(2l + 3)\\(2l-l)\\y Vl > 



az(2/ + 3)! 



r=0 



In the generalized case, ZRP with angular momentum I generates also 21 complex poles of 
the S- matrix, which correspond the quasi-stationary states (resonances). The DTs (|2.9| ) with 
the parameters b\ results in the S"-matrix elements for dressed GZRP 



exp(2i?7 ( 



a, 



-ik 2l+1 fbi-ik? 



di + ik 2l+1 \bi + ik J 

We can use Darboux transformation in order to add (or remove) poles of the S-matrix. 



3 Multi-center problem 

The principal observation allows to built a zero-range potential eigen function in the multi- 
center problem. Let us consider n ZRPs at the points r« and interaction The wavefunction 
^(r) can be expressed in terms of the (outgoing- wave) Green function, defined by the equation 

(-V 2 + 2v {1) - k 2 ) g+(r, r') = 4vr5(r - r') 

The second (ingoing- wave) Green function is defined by g~(r, r') = g + *(r, r'). The partial waves 
ty\(r, k), defined by 

*A(r, k) r -j^ (Ax(n)e tkr+ ^ + A x (-n) e - tkr -^) , (3.11) 

for multi-centered target which includes n ZRPs and interaction t>W can be expressed as a 
superposition of the Green functions 

n 

* A (r, k) = (2ik)~ l £ c Xl (e^g + (r, n) + e^^r, r 4 )) (3.12) 
i=i 
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in which rj\ are phase shifts, A\(n) denote S- matrix orthonormal eigenfunctions, c A j are real 
numbers. They naturally generalize the spherical partial waves (kr)~ 1 ipi(r, k)Yt m (n) for a non- 
spherical potential ||. Expanding the partial waves ty\(r,k) at the infinity we obtain the 
expressions for v4 A (n) 

A x (n)=jr CXl ¥ 1 \r l ,-kn). (3.13) 



i=l 



where *W(r, fen) r ~°° e ikn ' T + F^e ikr / r is wavefunction for potential t)W, and is scattering 
amplitude. Scattering amplitude for potential v^ 1 ' and n ZRP is given by 

F(n,n ) = F^(n,n ) + ^ £ ( c *» - l) A A (n)A A (n ). (3.14) 

A 



By the imposition of the boundary conditions ( |2.5|) the calculation of the partial waves is 
reduced to the solution of the following system 

X) c Ai (s( r i> r i) + c(r i ,r j )tan^ A ) = -c Ai (A; + 5^ + (a, + fox*) tanr? A ) , (3.15) 
in which we use 

g + (r, Ti) r ~ ri | — ^—r + + i(k + dki), (3.16) 

where s(rj,rj) = Img + (rj, r 3 -), c(rj,rj) = Reg + (r-j, rj). The tangents tan7/ A can be found from 
compatibility condition of this system. In simplest case v = 5oti = 5ki = 

_ cos(fcr^) _ sin(fcr^) _ _ 

'-'V i> jJ ' ^Mi "jj ? '"ij |m 

the system fl3.15j ) is reduced to the usual equations of ZRP theory ||. In order to construct 
functions s(Ti,Tj) and c(rj, r^) for dressed potential, i.e. for we need to write down Green 
function as single-center expansion over spherical harmonics. In the simplest case, when the 
prop function is (fo(r,ib) = sinh(fer), the Green function is given by 

g +( r y) _ ei * |r ~ r/| ^o\r,k)£\r',k) sin(Ar)e^ ^ < ^ 
|r-r'| krr' krr' 



where 



„m. 6coth(6r) - ik jkr ,m, ,. 6coth(6r) sin(Ax) - kcos(kr) 



Vk 2 + b 2 ru v ' ' \/FTF 

The function ■^q 1 '* (r, fc) has the following asymptotic at infinity 

ipo (r, k) r ^°° sin(fcr + 5), 5 = — arctan(/c/|6|). 

The integral cross section can be readily derived using the optical theorem: a = 47rfc~ 1 Im.F(no, n ) 
Thus, averaged integral cross section is given by 

a = ^sin(5) 2 + ^5:sin(r /A ) 2 . (3.18) 

A 
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4 X n and XY n structures 

For purpose of illustration we consider scattering problem for a dressed multi-center poten- 
tial. The multi-center scattering within the framework of the ZRP model was investigated 
by Demkov and Rudakov || (8 centers, cube), Drukarev and Yurova || (3 centers in line), 
Szmytkowski [|ll|] (4 centers, regular tetrahedron). 

4.1 Electron— X n scattering problem 

Let structure X n contains n identical scatterers, which involve only s waves. Denote position 
vector of the scatterer j by the rj. Suppose, for the sake of simplicity, the distance between any 
two scatterers i and j is r^- = R. There are three such structures in three-dimensional space - 



dome X2, regular trihedron X3, regular tetrahedron X4. The partial waves fl3.12|) and phases r]\ 
can be classified with respect to symmetry group representation, degeneracy being defined by 
the dimension of the representation ||. The structures X 2 , X 3 , X 4 belong to the A»7i> C31;, ?d 
point groups respectively. The equation (|3.15| ) leads to algebraical problem 



c\j(sm(kR) + ta.nr]\cos(kR)) = —c\i(kR + tan(r)\)aR), (4.19) 

3=1 

The phases can be readily found from compatibility condition of this system. To factorize the 
determinant we derive the expressions for the phases 

kR+(n- 1) sin(kR) 

tan ^i = S T\ tt^t, 4.20 

aR+{n-l) cos(kR) v ' 

kR - sm(kR) 



tanr M = p 7T^' A = 2,n. (4.21) 

aR — cos{kR) 



Thus, assuming n = 4 we obtain the phases of a regular tetrahedron JIT]]. The partial and 
integral cross sections can be expressed as 



47r tan(77 A ) 2 
k 2 1 + tan^) 1 



ax = TT l ; ^, 2 , a = a, + (n - l)a n . (4.22) 



4.2 Electron— YX n scattering problem 

The structures YX n can be used, for instance, to study a slow electron scattering by the 
polyatomic molecules like H 2 0, NH 3 , CH 4 , etc. Let ri,..,r n denote the position vectors of 
the scatterers X and r n+1 denotes the position vector of the scatterer Y. The scatterers X 
are situated in vertices of a regular structure X n , i.e. = R. Suppose, for the sake of 
simplicity, the distance between the scatterer Y and any scatterer X is rj( n +i) = D. Therefore 
the position of the scatterer Y perfectly fixed only if n = 4 (geometric center of the tetrahedron, 
R = 2y / f J D). The partial waves are given by the equation ( |3.12|) , where the summation should 
be performed over the i — l,n + 1. The constants c\i and phases can be derived analytically. 
Thus, we obtain the phases 

kR - sm(kR) 



t&nr]x = 71~d\ ' A = 3,n + 1, (4.23) 

aR — cos(kR) 



6 



60 
50 



C\J 

E 

cd° 40 
o 



30 



8 20 



10 




o 
o 

(J 

o 
o 


°°°°Oo 

°°o c 


cP 





4 6 
Impact Energy (eV) 



10 



Figure 1: Integral cross sections for electron— SiH4 
scattering. ZRP model: solid line; Experimental re- 
sults of Wan et al.[Q: open circles. 



where a is inverse scattering length of the scatterers X. The 
equation 



tan 771,2 obey the quadratic 



(k + f3x) x + 



kR+(n- 1) sm(kR) \ R(sm(kD) + x cos(kD)) 

; 7T T, „x = n——, ; : 7— - 



aR+(n- 1) cos(AlR) J D 2 (aR + (n - 1) cos(kR)) ' 



(4.24) 



where f3 indicates inverse scattering length of the scatterer Y. In the limiting case when the 
distance between X and Y scatterers is very large, i.e. D — > oo, the expression for tan 7/1 
becomes ( |4.20|) and tan7/ 2 — > —k/(3. This situation corresponds to independent scattering on 
the structure X n and scatterer Y. The tangent of 7/1 also reduces to tan rfi for structure X n in 
the limit a® co. 



4.3 Dressing of the structures X n and discussion 

The integral cross sections for SiH 4 (closed-shell ground-state 1 A\) are plotted in fig. [I] for a 
number of values of a,(3,R,D, which are regarded as constant in the range of interest. The 
present calculations were carried out with a = 0.33 clq , (3 = 0.41 Oq , R = 4.51 ao, D = 2.76 ao- 
Our calculations were made within the framework of the ZRP model and hence they are not 
expected to be correct for low impact energies (i.e for energies ~ 1 eV) where polarization 
effects (SiH 4 spherical polarizability ~ 30.4 a^) are known to be important. Because induced 
polarization potentials are always attractive, including polarization algebraically increases the 
computed phases (771,3 < at E ~ 0) and therefore decreases the integral cross section. Thus, 
fig. [I] clearly shows that polarization effects reproduce the deep minimum (Ramsauer-Townsend 
minimum near 0.35 eV) seen in the experimental data |12| and numerical calculations [|l0fl , 
which are incomparably smaller in the our calculations. The ICS can be corrected by the 
ZRPs dressing. Fig. ^| show ICS for dressed YX4 structure. In our calculation, the parameters 
a = 0.35, (3 = 0.38 were used. In the higher collision energies our cross section (fig. |l|) differs 
from the experimental results in size but coincides in shape. It is known that role of the higher 
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Figure 2: Integral cross sections for electron— S1H4 scat- 
tering on dressed ZRPs 

partial waves is important at the higher collision energies. Thus, taking into consideration the 
p waves for the scatterers X and Y we can add in partial waves and improve the agreement 
both in size at the higher energies and in the position of shape resonance. 

5 Conclusion 

We demonstrate the posibilities of DT in multi-center scattering problem. Thus, these trans- 
formations allow to correct the ZRP model at low energies. In the limiting case 6^0 DT 
induces the long-range forces ~ r~ 3 depending on angular momentum and leads to singularity 
in the cross section at zero energy. 
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